Constraints on a Cardassian model from SNIa data — revisited 



Wlodzimierz GodlowskQ and Marek Szydlowsk{j] 

Astronomical Observatory, Jagiellonian University, Orla 171, 30-244 Krakow, Poland 

Adam KrawiecH 

Institute of Public Affairs, Jagiellonian University, Rynek Glowny 8, 31-042 Krakow, Poland 

We discuss some observational constraints, resulting from SN la observations, imposed on the 
behavior of the original flat Cardassian model, and its extension with the curvature term included. 
We test the models using the Perlmutter SN la data as well as the new Knop and Tonry samples. 
We estimate the Cardassian model parameters using the best-fitting procedure and the likelihood 
method. In the fitting procedure we use density variables for matter, Cardassian fluid and curvature, 
and include the errors in redshift measurement. For the Perlmutter sample in the non-flat Cardassian 
model we obtain the high or normal density universe (fi m ,o > 0.3), while for the fiat Cardassian 
model we have the high density universe. For sample A in the high density universe we also find 
the negative values of estimates of n which can be interpreted as the phantom fluid effect. For 
the likelihood method we get that a nearly flat universe is preferred. We show that, if we assume 
O , that the matter density is 0.3, then n m in the flat Cardassian model, which corresponds to the 

Perlmutter model with the cosmological constant. Testing with the Knop and Tonry SN la samples 
' show no significant differences in results. 

o 

^ ! I. INTRODUCTION 

. Freese and Lewis |1| have recently proposed an alternative to the cosmological constant model explaining the 
currently accelerating Universe. In this model the standard Friedmann-Robertson- Walker (FRW) equation is modified 
| by the presence of an additional term p n , namely 
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H' = !- + Bp n , (1) 



' where H = a/a is the Hubble parameter, a is the scale factor, p is the energy density of matter and radiation, and B 



is a positive constant. For simplicity it is assumed that, the density parameter for radiation matter O r ,o = 0. This 



proposal seems to be attractive because the expansion of the universe is accelerated automatically due to the presence 
,H , of the additional term (if we put B = then the standard FRW equation is recovered) . 

Because the Cardassian model offers an alternative to the cosmological constant model,the agreement of this model 
with available observations of type la supernovae (SN la) was immediately verified 0,0. Some interesting results 
on observational constraints in the generalized Cardassian model were also obtained from the statistical analysis of 
gravitational lensed quasars 3 . In our analysis we use three SN la data sets compiled by Perlmutter et al. @ , Knop 
et al. [|| and Tonry et al. Q to test both the original Cardassian model and its extension with a curvature term. Our 
results are inconsistent with Zhu and Fujimoto 2j's prediction of the low density Universe. The noteworthy results 
from the joint analysis of SN la data and CMBR were obtained by Sen and Sen 0] who also confirmed the prediction 
of low density Universe, but their analysis was restricted to the case of n > 0. This case was also examined by Frith 
@ who used the Tonry SN la data. However, such a choice of an interval for this parameter is not physically justified 
and both negative and positive values are admissible. 

II. BASIC EQUATION 

We add a curvature term to the original Cardassian model. Then equation assumes the following form 

H* = P -+Bp n + ^. (2) 
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It is useful to rewrite equation in dimensionless variables fli t o, i = (m, Card, k), 

jj2 / a \ ^ / a \ 3n / \ 2 

772 = ^m,0 + ^Card,0 + ^fe.O : (3) 

Hfi \aoJ \a J \ a oJ 

where energy density for the dust is assumed (p = 0) (and hence p oc a -3 from a conservation equation), fl m< o is the 
matter density parameter; Slcard.o = 3Bp n /3Hq is the density parameter for the fictitious noninteracting fluid which 
mimics term p n and yields additional density in the model. The subscript zero denotes a present value of model 
parameters. Then 

p=-3H(p + p) (4) 
For a = ao (the present value of scale factor) we obtain the following constraint 

^m,0 + ^Card,0 + ^fc,0 = ^ W 

For the fictitious Cardassian fluid we obtain p — {n — l)p if p oc a~ 3 . In this interpretation the negative value of n is 
equivalent to the existence of phantom matter with supernegative presssure. 

Let us note that in the analogous equation describing FRW dynamics with dust and additional fictitious fluid with 
the equation of state p = (n — l)p presented in Sen and Sen Q only case n > was considered. 



III. MAGNITUDE-REDSHIFT RELATION 



It is well known that various cosmic distance measures, for example the luminosity distance, depend sensitively 
both on the spatial geometry (curvature) and dynamics. Therefore, luminosity depends on the present densities of 
the different components of matter content and its equation of state. For this reason, the magnitude-redshift relation 
for distant objects is proposed as a potential test for cosmological models and plays an important role in determining 
cosmological parameters. 

Let us consider an observer located at r = at the moment t = to who receives light emitted at t = t± from a 
source of absolute luminosity L located at the radial distance r\ . It is known that the cosmological redshift z of the 
source is related to t\ and tg by the relation 1 + z = a(to)/a(ti). If the apparent luminosity / of the source measured 
by the observer is defined by 

then the luminosity distance of the source is 

d L = (1 + z)a Q r 1 . (7) 

For historical reasons, the observed and absolute luminosities are defined in terms of K-corrected observed magni- 
tudes and absolute magnitudes m and M, respectively (I = 10~ 2m / 5 • 2.52 • 10~ 5 erg cm -2 s~ 2 , L = 10~ 2M / 5 • 3.02 • 
10 35 ergs -2 ) 0. When written in terms of m and M, equation @ yields 

m(z, M, fi m ,o, ^Card : o, fife.o) = M + 5 log w [D L (z, f2 m , , £!card : o, fife.o)] (8) 

where 

M = M - 5 log 10 H + 25 (9) 

and 

~Dl(z, fi mj o, ^Card.O, ^/c,o) = H dL(z, fl m fi, ^Card,0, ^k,Q, H ) (10) 

is the dimensionless luminosity distance while dh is in Mpc. 

The standard analysis yields the following relationship for the dimensionless luminosity distance 

2? L ((«,n m ,o,nc»«i,o,nfc,o) = ^ / [(i-a„,o-^ca rd ,o)(i + ^) 2 

+ Q m , (l + z'f + fi C ard,0(l + z'? n r 1/2 dz' I , (11) 
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TABLE I: Results of the statistical analysis of the Cardassian model obtained both for Perlmutter sample A and C from the 
best fit with minimum \ 2 (denoted with BF) and from the likelihood method (denoted with L). The same analysis was repeated 
with fixed f2 m ,o = 0.3. 
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where 



£(x) = sin a; with K, = —£lk,o when &k,0 < 

£(x) = x with JC = 1 when Q kfi = (12) 

£(:c) = sinha; with K, = £lk,o when il^.o > 



and the density parameter for hypothetical curvature fluid is 

k 



Thus, for given M., f2 m o; ^Card.O: ^fe,0) equations JHJ and (jllfl give the predicted value of m(z) at a given z. 
The goodness of fit is characterized by the parameter 

* = 2 ^2 +a 2 (13) 

where /j,q i is the measured value, /_tg i is the value calculated in the model described above, <r^ i is the measurement 
error, a^ z i is the dispersion in the distance modulus due to peculiar velocities of galaxies. 

We assume that supernovae measurements come with uncorrelated Gaussian errors and in this case the likelihood 
function C can be determined from a chi-squared statistic C oc exp (— x 2 /2) 



IV. STATISTICAL ANALYSIS WITH THE PERLMUTTER SAMPLE 

We have decided to test our model using the Perlmutter samples of supernovae . We estimate the value of M. 
(equation 10) separately for the full sample of 60 supernovae (sample A) and for the sample of 54 supernovae (sample 
C — in that sample we exclude 2 supernovae as outliers and 2 as likely reddened ones from the sample of 42 high 
redshift supernovae and 2 outliers from the sample of 18 low redshift supernovae). We obtain the value of M. = —3.39 
for sample A, and M. = —3.43 for sample C. To test the model we calculate the best fit with minimum x 2 as well as 
estimate the model parameters using the likelihood method ^3 • For both statistical methods we take the parameter 
ri in the interval [—3.33,2], f2 nii o m the interval [0,1], £lk,o in the interval [—1,1], while an interval for f2card,o is 
obtained from equation JSJ). 

First, we estimate best fitting parameters for the non-flat Cardassian model and obtain ilk.o = —1.00, f2 m ,o = 0.83, 
n = —0.67, £!card,o = 1-17 for sample A, and Qk,o = —1-00, il mj o = 0.60, n — 0.50, f^card.o = 1-40 for sample C. This 
model is characterized with strong negative curvature, that is similar to the Perlmutter results. Let us note that such 
a high value of £lk,o is in contradiction with CMB observations, which prefer the flat universe. When we assume that 
^fc.o = then we obtain the best-fitted flat Cardassian model with O m o = 0.52, n = —0.93, f2card,o = 0.48 for sample 
A, and f2 m ,o = 0.46, n = —0.60, ficard.o = 0.54 for sample C. The detailed values are in Tables HI and ILT1 

For sample C we draw the best- fitted, non-flat (upper- middle curve), flat (lower-middle curve) Cardassian models as 
well as the best-fitted Perlmutter model with f2 m .o = 0.28, JIa.o = 0.72 (upper curve) against to the flat Einstein-de 
Sitter model (zero line) on the Hubble diagram (Fig. [TJ. One can observe that the difference between the best- 
fitted Cardassian model and the Einstein-de Sitter model with flA,o = assumes the largest value for z ~ 0.9 and 
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TABLE II: Results of the statistical analysis of the Cardassian flat model for Perlmutter sample A and sample C obtained both 
from the best fit with minimum \ 2 (denoted with BF) and from the likelihood method (denoted with L). In the case in which 
we marginalize over M we denote it with M. The same analysis was repeated with fixed f2 m ,o = 0.3. 
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TABLE III: Non-flat Cardassian model parameter values obtained from the minimization procedure carried out on the Perl- 
mutter samples. 



sample fl k ,o ficard,o 



m.O 



A -0.12t°i* 0.37t°f 0.50l° f 2 -0.13t?f 7 

c o.i7±g:SS o.33lg-» o.29±8:i| o.wtlit 



significantly decreases for higher redshifts. While the best-fitted flat Cardassian model and Perlmutter one have 
increasing differences to the flat Einstein-de Sitter model for higher redshifts, the differences between the best-fitted 
flat Cardassian model and Perlmutter model increase for higher redshifts. It gives us a possibility to discriminate 
between the Perlmutter model and Cardassian models if z ~ 1 supernovae data will be available. It is very important 
because for the present data the Cardassian models are only marginally better than the Perlmutter model. 

However, knowing the best-fit values alone has no enough scientific relevance, if confidence levels for parameter 
intervals are not presented, too. Therefore, we carry out the model parameters estimation using the minimization 
procedure, based on the likelihood method. On the confidence level 68.3% we obtain parameter values for samples A 
and C (Table [Trill. 

It should be noted that values obtained in both methods are different but the minimization procedure seem to be 
more adequate for analyzing our problem. The density distributions /(rJfc,o) and fiji) in the Cardassian model are 
presented on Fig. [3 and |3 respectively. These figures show that the preferred model of the Universe is the nearly 
flat one, that is in agreement with CMBR data. For sample C the most probable value of f2 m ,o is 0.29 that is in 
agreement with the present CMBR and extragalactic data [jj], [T^] • Therefore, the detailed analysis of the flat case of 
the Cardassian model is carried out. 

In the Fig. 2| we presented confidence levels on the plane (^ m ,0j n ) minimized over M. for the flat Cardassian model. 
Figure Q] shows the preferred values of parameters f2 mi o and n. We find that the expected value of f2 m .o increases 
when n decreases. The similar result was already obtained by Sen and Sen |3J, for n > 0. It should be pointed out 
that these authors argued that f2 mj o could be less than 0.3 from the confidence level obtained from the joint analysis 
of supernovae data and the positions of the CMBR peaks. In the Cardassian model the structure formation would 
be significantly different from, e.g., the classical model with the cosmological constant. Consequently, in order to 
calculate the consequence for the CMB the simple angular rescaling does not suffice. 

The detailed results of our analysis for the flat model are summarized in Table ITU The best fit procedure suggests 
that n should be negative and consequently f2 m .o is greater than 0.3. While we find the parameter n is negative for 
minimum value of % 2 statistic, the confidence levels allows the positive values of n. For the maximum likelihood method 
on sample A we obtain that fl m ,o — 0.60t^i5 ano ^ n = ~ ^-^-{.27 011 the confidence level 68.3% for M. = —3.39; while 
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FIG. 1: Residuals (in mag) between the Einstein-de Sitter model and four cases: the Einstein-de Sitter itself (zero line), 
the Perlmutter flat model (upper curve), the best-fitted flat Cardassian model (upper-middle curve), Qk.o = 0, f2 m ,o = 0.46, 
n = —0.60, ficard.o = 0.54 and the best-fitted (non-flat) Cardassian model (lower-middle curve) flk,o ~ —1.00, £7 m ,o = 0.60, 
n = 0.50, ficard.o = 1-40. 

^m,o = 0.54to'ii and n = — 1.20^5 '|g on the confidence level 68.3% when we marginalize over M. 

In turn for sample C we obtain that Q^q = 0.58t° °il and n = O- 10 ^? on tne confidence level 68.3% for M = -3.43 
(Fig- El and [fj), while O m o = 0.52+q'i3 and n = 0. 17^2 04 011 the confidence level 68.3% when we marginalize over M. 

For the flat model with f2 m ,o = 0.3 we obtain for sample A n = —0.03 with a(n) = 0.13 both for M. — —3.39 and 
with <j{n) = 0.33 when we marginalize over M. In turn for sample C we obtain that n = 0.03 with a(n) — 0.12 both 
for M = —3.43 (Fig.0 and <r(n) = 0.23 when we marginalize over M. 

V. STATISTICAL ANALYSIS WITH THE KNOP SAMPLE 

Because the Perlmutter sample was completed four years ago, it would be interesting to use newer supernovae 
observations. Lately Knop et al. Q have reegzamined the Permutter sample with host-galaxy extinction correctly 
applied. They chose from the Perlmutter sample these supernovae which were the more securely spectrally identified 
as type la and have reasonable color measurements. They also included eleven new high redshift supernovae and a 
well known sample with low redshift supernowae. 

We have also decided to test the Cardassian model using this new sample of supernovae. They distinguished a 
few subsets of supernovae from this sample. We consider two of them. The first is a subset of 58 supernovae with 
extinction correction (Knop subsample 6; hereafter K6) and the second one of 54 supernovae with low extinction (Knop 
subsample 3; hereafter K3). Sample C and K3 are similarly constructed because both contain only low extinction 
supernovae. 

It should be pointed out that in contrast to Perlmutter et al. [5| who included errors in measurement of redshift z 
Knop et al. |6| took only uncertainties in the redshift due to peculiar velocities. It is the reason that for security we 
separately repeated our analysis including errors in measurement of redshift z (subsamples K6z and K3z) . The errors 
was taken from Perlmutter et al. In the case when errors were not available we assume a(z) = 0.001. As one can 
see above, this change has only marginal influence on the results. 

At first we estimate the value of M. (equation ©) separately for both sample K6 and K3. We obtain the value of 
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FIG. 2: The density distribution for O fc>0 in the Cardassian model (sample C). We obtain that O fei o = 0.171^66 on tne confidence 
level 68.3% (the inner dash lines). Both positive and negative values of Slk,o are formally possible. 
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FIG. 3: The density distribution for n in the Cardassian model (sample C). We obtain that n = O.IOlt^g on the confidence 
level 68.3% (the inner dash lines). Both positive and negative values of n are formally possible. 
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FIG. 4: Confidence levels on the plane (fi m ,o,n) minimized over for the flat model (£lk,o = 0), and with ficard.o = 1 — fim,o- 
The figure shows of the preferred value of £7 m ,o and n. 
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FIG. 5: The density distribution for Q m ,o in the Cardassian flat model. We obtain that S7 m ,o = 0.581q;j3 on the confidence 
level 68.3% (the inner dash lines). Additionally the confidence level 95.4% is marked (the outer dash lines). 
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FIG. 6: The density distribution for n in the Cardassian flat model. We obtain that n = O.lOl?;^ on the confidence level 68.3% 
(the inner dash lines). Additionally the confidence level 95.4% is marked (the outer dash lines). Both positive and negative 
values of n are formally possible. 
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FIG. 7: The density distribution for n in the Cardassian flat model with fi m ,o = 0.3. We obtain that n — 0.03 with a(n) = 0.12. 
Both positive and negative values of n are formally possible. 
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TABLE IV: Results of the statistical analysis of the Cardassian model obtained both for the Knop samples from the best fit 
with minimum \ 2 (denoted with BF) and from the likelihood method (denoted with L). The same analysis was repeated with 
fixed fi m ,o = 0.3. 
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jM = —3.53 for sample K6, and M. = —3.48 for sample K3. It is in agreement with Knop et al. 0. 

First, we estimate best fitting parameters for non-flat Cardassian model and obtain fi^ q = —0.24, f2 m = 0.65, 
n = —1.10, f2card,o = 0-59 for sample K6, and fi fe-0 = —0.36, O m ,o = 0.49, n — —0.20, Qcard.o = 0.87 for sample K3. 
When we assume that Clkfi — (flat universe) then we obtain the best-fitted flat Cardassian model with £l m .o = 0.53, 
n = —1.50, f^card.o = 0-47 for sample K6, and f2 mj o = 0.41, n = —0.63, ^card.o = 0-59 for sample K3. The detailed 
values are in Tables ITVI and M 

We also carry out the model parameters estimation using the minimization procedure, based on the likelihood 
method. On the confidence level 68.3% we obtain parameter values for samples K6 and K3 (Table lvi*|) . 

It should be noted that values obtained in both methods are different but differences are much smaller then in the 
case of original Perlmutter sample. 

The detailed results of our analysis for the flat model are summarized in Table [V] The best fit procedure suggests 
that n should be negative and consequently f2 mj o is greater than 0.3. While we find the parameter n is negative for 
minimum value of x 2 statistic, the confidence levels allows the positive values of n. For the maximum likelihood method 
on sample K6 we obtain that fl m ,o = 0.44+Q'Jg and n = — 1.13^^3 on the confidence level 68.3% for M. = —3.53; 
while f2 m .o = 0.48loog an d n = — 3.33 +2 ° on the confidence level 68.3% when we marginalize over M. 

In turn for sample K3 we obtain that fim.o = 0.50+q 15 an( i n = — 0.37^°'qq on the confidence level 68.3% for 
M = —3.48, while Q ro> o = 0.48to'i3 an d n = ~ O^ti'll on the confidence level 68.3% when we marginalize over M. 

For the flat model with Sl nii o = 0.3 we obtain for sample K6 n = —0.17 with <j(n) = 0.16 for Ai = —3.39 and 
71 = —0.13 with cr(n) — 0.23 when we marginalize over A4. In turn for sample K3 we obtain that n = —0.20 with 
cr(ra) = 0.13 for Ai = —3.48 and n = —0.20 u{n) = 0.17 when we marginalize over Ai. 

VI. STATISTICAL ANALYSIS WITH THE TONRY SAMPLE 

Another sample was presented by Tonry et al. 7] who collected a large number of supernovae published by different 
authors and added eight new high rcdshift SN la. This sample of 230 SNe la was recalibrated with consistent zero 
point. Whenever it was possible the extinctions estimates and distance fitting were recomputed. However, none of 
the methods was able to apply to all supernovae (for details see Tab. 8 0)- 

Despite of this problem, the analysis of the Cardassian model using this sample of supernovae could be interesting. 
We decide to analyse four subsamples. First, the full Tonry sample of 230 SNe la (hereafter sample Ta) is considered. 
The sample of 197 SNe la (hereafter sample Tb) consists of low extinction supernovae only (median V band extinction 
Ay < 0.5). Because the Tonry sample has a lot of outliers especially in low redshift, the sample of 195 SN la is such 



10 



TABLE V: Results of the statistical analysis of the Cardassian flat model for the Knop samples obtained both from the best fit 
with minimum \ 2 (denoted with BF) and from the likelihood method (denoted with L). In the case in which we marginalize 
over M we denote it with M. The same analysis was repeated with fixed Q m ,o = 0.3. 



sample 


- " ^tLI (J. , U 


n o 


n 


M 


x 2 


method 


K6 


0.47 


0.53 


-1.50 


-3.53 


54.6 


BF 




0.44 


0.56 


-1.13 


-3.53 




L 




0.46 


0.54 


-3.33 


-3.60 


53.5 


M, BF 




0.48 


0.52 


-3.33 


-3.55 




M, L 




0.70 


0.30 


-0.17 


-3.53 


55.7 


BF 




0.70 


0.30 


-0.17 


-3.53 




L 




0.70 


0.30 


-0.10 


-3.52 


55.6 


M, BF 




0.70 


0.30 


-0.13 


-3.53 




M, L 


K6z 


0.45 


0.55 


-1.87 


-3.53 


43.6 


BF 




0.43 


0.57 


-1.57 


-3.53 




L 




0.46 


0.54 


-3.33 


-3.60 


42.7 


M, BF 




0.48 


0.52 


-3.33 


-3.54 




M, L 




0.70 


0.30 


-0.17 


-3.53 


44.9 


BF 




0.70 


0.30 


-0.17 


-3.53 




L 




0.70 


0.30 


-0.07 


-3.51 


44.8 


M, BF 




0.70 


0.30 


-0.10 


-3.51 




M, L 


K3 


0.59 


0.41 


-0.63 


-3.48 


60.3 


BF 




0.50 


0.50 


-0.37 


-3.48 




L 




0.58 


0.42 


-0.77 


-3.49 


60.3 


M, BF 




0.52 


0.48 


-0.40 


-3.49 




M, L 




0.70 


0.30 


-0.20 


-3.48 


60.4 


BF 




0.70 


0.30 


-0.20 


-3.48 




L 




0.70 


0.30 


-0.13 


-3.47 


60.4 


M, BF 




0.70 


0.30 


-0.17 


-3.48 




M, L 



K3z 0.60 0.40 -0.57 -3.48 52.0 BF 

0.49 0.51 -0.30 -3.48 L 

0.58 0.42 -0.83 -3.50 52.0 M, BF 

0.52 0.48 -0.40 -3.50 M, L 

0.70 0.30 -0.17 -3.48 52.2 BF 

0.70 0.30 -0.17 -3.48 L 

0.70 0.30 -0.17 -3.48 52.2 M, BF 

0.70 0.30 -0.17 -3.48 M, L 



TABLE VI: Model parameter values obtained from the minimization procedure carried out on the Knop samples. 



sample ft k . ricard.o ^m,o n_ 

K6 -0.02±° ;« 0.39±S^ 0.45l° j? -0.431°;™ 

K6z -0.15l°i« 0.431°-^ 0.52l^o -0.50t°/g 

K3 0.29j£» 0.34t°-^ 0.27±° ^ -0.031?!? 

K3z 0.30l° ^ 0.34l° « 0.25!°,^ -O.OOl?'* 



that all low redshift (z < 0.01) supernovae are excluded (hereafter sample Tc). In the sample of 172 SN la all 
supernovae with low redshift and and high extinction arc omitted (hereafter sample Td). 

Tonry et al. [3 presented redshift and luminosity distance observations for their sample of supernovae. Therefore, 
equations (JSJ and Jj|J should be modified }13|: 

TO - M = 5 log 10 (Z> L ) Tonry - 5 log 10 65 + 25 (14) 
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TABLE VII: Results of the statistical analysis of the Cardassian model obtained for the Tonry samples both from the best fit 
with minimum \ 2 (denoted with BF) and from the likelihood method (denoted with L). The same analysis was repeated with 
fixed f2 m ,o = 0.3. 



sample 


^Card 


^fc 





n 


M 


x 2 


method 


Ta 


1.52 


-1.00 


0.48 


0.27 


15.935 


252.1 


BF 




0.54 


-0.64 


0.49 


0.10 


15.935 




L 




1.68 


-0.98 


0.30 


0.37 


15.935 


252.5 


BF 




0.28 


0.42 


0.30 


0.33 


15.935 




L 


Tb 


1.50 


-1.00 


0.50 


0.23 


15.935 


185.3 


BF 




0.64 


-0.78 


0.49 


0.07 


15.935 




L 




1.70 


-1.00 


0.30 


0.37 


15.935 


185.7 


BF 




0.31 


0.39 


0.30 


0.33 


15.935 




L 


TV 


1 47 


- l.UU 


0.48 


0.23 


15.935 


ZULU 






0.54 


-0.64 


0.49 


0.10 


15.935 




L 




1.68 


-0.98 


0.53 


0.37 


15.935 


201.3 


BF 




0.28 


0.42 


0.49 


0.33 


15.935 




L 


Td 


1.50 


-1.00 


0.50 


0.23 


15.935 


164.9 


BF 




0.63 


-0.78 


0.49 0.17 15.935 




L 




1.70 


-1.00 


0.30 


0.37 


15.935 


165.4 


BF 




0.31 


0.39 


0.30 


0.33 


15.935 




L 



and 

M = -51og 10 iJo + 25. (15) 

For H = 65 km s -1 Mpc -1 wc obtain M = 15.935. 

First, we estimate best fitting parameters for the non-flat Cardassian model and obtain f^.o = — 1, ^ m .o = 0.48, 
n = —0.27, ricard.o = 1-52 for sample Ta. When we assume that fi^o = (fiat universe) then for the sample Ta we 
obtain the best-fitted flat Cardassian model with r2 mj o = 0.46, n = —1.50, Slcard.o = 0.54. For samples Tb, Tc, and 
Td results differs only marginally. The detailed values are presented in Tables IVIII and IVIIII 

As in previous sections, we also carry out the model parameters estimation using the minimization procedure based 
on the likelihood method. On the confidence level 68.3% we presented obtained parameter values for Tonry samples 
in the Table E| 

The detailed results of our analysis for the flat model, based on the likelihood method, are summarized in Table IVlTTl 
For sample Ta we obtain that Q mi0 = 0.51±£?2 and n = -°- 40 tcK63 on the confidence level 68.3% for M = 15.935; 
while f2 m ,o = 0.53ig'og and n = — 1.87ii"oo on ^ ne confidence level 68.3% when we marginalize over M.. 

In turn for sample Td we obtain that il mj o = 0.491q'°2 and n = — 0.43tp'g3 011 the confidence level 68.3% for 
M = 15.935, while f2 m ,o = 0-5liolo6 and n = — l'2Q-iI)6 on tne confidence level 68.3% when we marginalize over M. 

For the flat model with fi m ,o = 0.3 we obtain for subsample Ta n = 0.03 with <j{n) = 0.10 for M. = 15.935 and 
n = —0.10 with o(n) — 0.15 when we marginalize over M.. For other subsamples results are very similar. 

All these results suggest that n is negative and consequently r2 m ,o is greater than 0.3. 

Generally we find that results obtained with the Tonry sample are the similar as obtained with the Perlmutter 
sample apart from the non flat case where values of n parameter are positive, because either positive and negative 
values are in the Perlmutter samples. 

VII. DISCUSSION 

In the present paper we discussed the problem of universe acceleration in the Cardassian model. Our results are 
different from those obtained by Zhu and Fujimoto who suggested the universe with very low matter density. Our 
results indicate the high or normal (Sl m ,o ~ 0.3) matter density. This difference may come from not including the 
errors in redshifts and using a different variable by them. In our fitting procedure we use the more natural variable 
^Card.o (or ^ m ,o) instead of their parameter z cq (z cq : | = Bp n , p = poa~ 3 , 1 + z = a -1 ). Two parameters fl m .o and 
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TABLE VIII: Results of the statistical analysis of the Cardassian flat model for the Tonry samples obtained both from the best 
fit with minimum \ 2 (denoted with BF) and from the likelihood method (denoted with L). In the case in which we marginalize 
over M. we denote it with M. The same analysis was repeated with fixed f2 m ,o = 0.3. 



sample 


V_/ ell CI , U 





n 


M 


x 2 

A, 


method 


Ta 


0.54 


0.46 


-0.50 


15.935 


253.6 


BF 




0.49 


0.51 


-0.40 


15.935 




L 




0.48 


0.52 


-2.00 


15.855 


247.7 


M, BF 




0.47 


0.53 


-1.87 


15.855 




M, L 




0.70 


0.30 


0.03 


15.935 


254.7 


BF 




0.70 


0.30 


0.03 


15.935 




L 




0.70 


0.30 


-0.10 


15.895 


252.1 


M, BF 




0.70 


0.30 


-0.10 


15.895 




M, L 


Tb 


0.56 


0.44 


-0.50 


15.935 


187.1 


BF 




0.51 


0.49 


-0.40 


15.935 




L 




0.51 


0.49 


-1.47 


15.875 


183.4 


M, BF 




0.49 


0.51 


-1.40 


15.875 




M, L 




0.70 


0.30 


-0.03 


15.935 


188.1 


BF 




0.70 


0.30 


-0.03 


15.935 




L 




0.70 


0.30 


-0.13 


15.905 


186.5 


M, BF 




0.70 


0.30 


-0.17 


15.905 




M, L 



Tc 


0.54 


0.46 


-0.50 


15.935 


202.5 


BF 




0.49 


0.51 


-0.40 


15.935 




L 




0.49 


0.51 


-1.53 


15.875 


198.7 


M, BF 




0.47 


0.53 


-1.47 


15.875 




M, L 




0.70 


0.30 


0.03 


15.935 


203.6 


BF 




0.70 


0.30 


0.03 


15.935 




L 




0.70 


0.30 


-0.07 


15.905 


202.2 


M, BF 




0.70 


0.30 


-0.07 


15.905 




M, L 


Td 


0.56 


0.44 


-0.50 


15.935 


166.7 


BF 




0.51 


0.49 


-0.43 


15.935 




L 




0.52 


0.48 


-1.23 


15.885 


164.1 


M, BF 




0.49 


0.51 


-1.20 


15.885 




M, L 




0.70 


0.30 


-0.03 


15.935 


167.7 


BF 




0.70 


0.30 


-0.03 


15.935 




L 




0.70 


0.30 


-0.13 


15.905 


166.8 


M, BF 




0.70 


0.30 


-0.13 


15.905 




M, L 



TABLE IX: Model parameter values obtained from the minimization procedure carried out on the Tonry sample. To obtain 
errors of the parameter flk,o we enlarge an estimation interval of flk,o to [—2, 2]. 



sample Q k . fl C ard,o ^m,o n_ 

Ta -0.641";™ 0.54l° j 7 0.49l° ;£ 0.10l°;f 

Tb -0.781°;™ 0.64l° j° 0.491°;^ 0.07l°; 3 6 ? 



Tc -0.64+-- 0.54-J;- 0.491-- 0.10+^0 

Td -0.781°'" 0.631" ^ 0.491°-^ 0.17l° f 7 



z cq are related by the formula 

^m,0 = (l + (l + ^c q ) 3(1 ^ ) )- 1 . (16) 
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FIG. 8: Relation between Qm,o and z eq for different n. 



In Fig.|Hlrelation (jlfijl for different n is presented. We can observe that taking a constant step in z eq it leads to varying 
steps in il m o- I n our method we have 

^m,0 = 1 — ^Card,0 — ^fc,0 (17) 

and we take the constant step of 0.01 in flm,o- Fot the flat Cardassian model both Zhu and Fujimoto 0] and we 
obtained negative values of n with the best fit method. 

Assuming the curvature term flk.o in the Cardassian model we obtained its statistical estimation to be close to 
zero. In this natural way we got that a nearly flat universe is preferred. Moreover, in this case for Pcrlmuttcr sample 
C the preferred value of fl m .o is 0.29 that is in agreement with CMBR and extragalactic data. While for sample A 
we got that a universe with higher matter content is preferred. 

Let us note that Perlmutter et al. Riess et al. [13 obtained the high negative value of £lk,o as the best fit, 
although zero value of flk,o is statistically admissible. To find the curvature they additionally used the data from 
CMBR and extragalactic astronomy. But in the Cardassian model in natural way the nearly flat universe is favored. 

If the flat universe is assumed, in the analysis of the Cardassian universe with sample C we got low but positive 
values of n and matter density higher than 0.3. For all other samples we obtained negative values of n (it would also 
indicate the phamtom fluid with super negative pressure). In turn, when £l m< o is fixed to be 0.3, then n ~ is favored. 

In general, the results obtained using the Perlmutter and Knop samples are not significantly different. The main 
advantage of using the Knop sample is the lower errors of estimated parameters in the model. The Tonry sample 
also gave the results similar to the Perlmutter sample. It can be noted that the Tonry sample does not improve the 
estimation errors compared with the Perlmutter sample. 

From the statistical analysis we found that the Perlmutter and Cardassian models are indistinguishable. Hence, 
using the philosophy of Occam's razor the Cardassian model should be rejected if we base on current SN la data. 
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However, the Cardassian model has its own advantages and fits well to the present observational data. And till we 
have no more precise data we should wait before we reject this model. 
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